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Summary 


This report reviews initial progress on a new program of investigation 
of scintillation and coherence effects for a laser-illuminated , noncoopera¬ 
tive target as viewed through atmospheric turbulence. The application of 
the effort is in the prediction of turbulence effects on the operation of 
coherent optical adaptive transmitter (COAT) systems. 

Significant analytical progress is reviewed, including results for the 
mutual coherence function, variance and covariance of irradiance, spectra, 
and statistics for a coherently illuminated diffuse target. The approach 
utilizes the extended Huygens-Fresnel principle, and includes turbulence 
effects in both the path from the transmitter to the target and that back 
to the receiver. It is found that there are three pertinent covariance 
sc • “s and six possible parameter realms, and that the normalized variance 
will t>e unity except in those cases when the target spot is sufficiently 
omall as to constitute a quasi-point-source. This parameter-realm view¬ 
point is further explored in relation to a variety of possible sources 
operating through turbulence. The analysis is partially extended to a 
more complex target and to the incoherent case, and future analytical tasks 
and applications to real adaptive systems are outlined. 

The establishment of an experimental field facility is also described, 
which will be capable of measuring all pertinent quantities at both visible 
and middle-infrared wavelengths./ Preliminary experimental results are 
presented. 

The results of this work should aid materially in the understanding of 
the dynamic behavior of adaptive laser sources operating in the presence of 
significant atmospheric turbulence. 
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I. Introduction 

This report reviews a new and ongoing investigation of scintillation 
effects arising from the illumination of noncooperative targets with laser 
radiation in the presence of atmospheric turbulence. The application of 
interest is an understanding of the performance of adaptive optical transmit¬ 
ter systems in the presence of target and turbulence-induced speckle and 
other scintillation structure. 

During the initial period of research, a theoretical understanding of 
turbulence effects was achieved for diffuse targets with coherent and 
incoherent (highly-multimode) illumination'. Progress was also made on 
the problem of target structure (glints) ; and on scintillations from gen¬ 
eralized sources, including parameter realms and covariance scales. These 
theoretical results are discussed in detail in Section II. 

During the same period, an experimental program was developed for 
visible and infrared (3.5-3.8 micron) wavelengths. Preliminary field 
experiments were conducted at 6328^ and a well-instrumented field facility 
established at 4880&. Concurrently, a pulsed infrared laser and receiver 
system was designed and fabrication initiated. The experimental efforts 
are discussed 'n Section III. 

II. Theoretical Description 

Significant progress has been made in the physical and analytical 
understanding of turbulence effects cn the dynamics of target-reflected 
radiation. The important quantities, relating ultimately to the performance 
of an adaptive transmitter system, are the variance of irradiance (o^), 
covariance of irradiance [C^(p)], mutual coherence function [T(p)], 
probability distribution function of irradiance, and spatial and temporal 
power spectrum. These quantities are of interest in the plane of the active 
laser transceiver, and include turbulence effects on both the illuminating 
radiation from the transmitter to the target and on the scattered radiation 
over the return path (Fig. 1). 

The treatment to be given below utilizes primarily the extended Huygens- 
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1 2 

Fresnel principle as applied to a turbulent path. ’ Until otherwise 
stated, the laser is assumed to be a coherent (TEM ) source, collimated 

X»U 

or focused, with a perfectly diffuse target. In all cases, we attempt to 
show clearly the assumptions and approximations that are made and to dis¬ 
cuss their implications. An important feature is the definition of dis¬ 
tinct (asymptotic) parameter conditions applying to any given configuration; 
there are generally six such conditions, representing the possible permuta¬ 
tions of inequalities between the three pertinent parameters: Fresnel zone 
% 

size (L/k) , coherence radius (o ) and speckle size in the absence of 

o 

turbulence. Each such condition will in general carry a distinct physical 
and analytical interpretation. It will be seen that cases of strong scin¬ 
tillations ("saturation" or multiple scattering) are included in these con¬ 
ditions, so that the treatment is general. 

We first consider (Sec. II-A) the most common situation, i.e., that in 
which the primary effect of the turbulence on the reflected radiation 
arises through the perturbation of the phase term in the associated Green's 
function. The implications of this assumption in terms of the field sta¬ 
tistics are also explained. A review of some of the developments of this 

3 

section appeared in the preceding technical report but without a full 
interpretation of the assumptions. 

We then generalize the development (Sec. II-B) to include the effects 
of the amplitude perturbation term in the Green's function. The implica¬ 
tions are again explored in detail. In Sec. II-C, we treat the first-order 
statistics of irradiance for a target containing one or more glints. 
Incoherent illumination is then discussed in Sec. II-D. 

In Sec II-E, we relate the current development to an alternative, 
earlier treatment of a related problem and show that the two approaches 
are basically complementary. We also consider the possible parameter realms 
for a variety of source types, with attendant predictions of strengths and 
spatial scales of scintillations. 

Finally, in II-F, we outline further topics and extensions of interest 
and related efforts to be undertaken. 

1. R.F.Lutomirski and H.T.Yura,"Propagation of a Finite Optical Beam in an 
Inhomogeneous Medium", Applied Optics.H), 1652, July 1971. 

2. H.Yura,"Mutual Coherence Function of a Finite Cross Section Optical Beam 
Propagating in a Turbulent Medium",Applied Optics, 1399, June 1972. 

3. J.R.Kerr, et al, "Propagation of Multiwavelength Laser Radiation through 
Atmospheric Turbulence", RADC Technical Report, August 1975. 
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II.A.Basic Irradiance Statistics and Mutual Coherence Function 


Previous work on speckle statistics has primarily been concentrated 

on the nature and statistics of the target surface, propagation of the 

speckle field without turbulence, and effects of speckle on image quality. 

Speckle propagation through turbulence has been considered over a vertical 

A 

path for the purposes of speckle interferometry. 

In the present section an analysis is given of the first and second 
order statistics of the received intensity (irradiance) after scattering 
from a diffuse target. The treatment is based on the extended Huygens- 
Fresnel formulation and includes the effects of the turbulent atmosphere 
on the laser beam as it propagates to the target and on the speckle as it 
propagates back to the receiver, Formulations aru- given for both the focused 
and collimated cases. The analysis also includes the mutual coherence func- 
tion (MCF). 

The source, target, and receiver configuration is shown in Figure 1. 

The present analysis is confined to the case of a TEM Q0 laser illuminator. 

The source and target are assumed to be much smaller than the path length 
(L), and the distance between the receiver and source is greater than the 
source size and much smaller than the path length. These geometric con¬ 
ditions confine the problem to small angles and ensure that the outgoing 
and returning radiation experience independent turbulence regions; the 
latter limitation is thought to be inessential owing to the diffuse target 
characteristics. 

1. Mean Irradiance at Receiver 

To find the mean irradiance, we need no assumptions other than 
that of a diffuse target. 

We write the source amplitude distribution as 


U (r) = U exp 
o o 


r 2 

2 a 2 


i kr 2 
' 2F 


( 1 ) 


and F are the characteristic beam radius and focal length respec- 

4. M. Elbaum, et al, "Laser Correlography* Transmission of High-Resolution 
Object Signatures through the Turbulent Atmosphere", Riverside Research 
Institute, Technical Report T-l/306-3-11, October 31, 1974. 



tively. The field at the target is written from the extended Huygens- 

1 2 

Fresnel principle ’ as 


U(PJ 


ke 

2 


ikL r 

k- /», 


(r) exp 


ik |p ~ r 1 2 + 

2L 


(P,r) 


dr 


(2) 


where describes the effects of the random medium on the propagation 
of a spherical wave. Combining Eqs. (1) and (2), we have 



o 


+ >,'•] (P,r) 


dr 


( 3 ) 


In particular, this applies to the special cases of a focused (L = F) 
or collimated (F -*■ °°) beam respectively. 

The field at the receiver is written by reapplying the Huygens- 
Fresnel principle to the field at the '.arget: 

ik [L + p£'| 

U (P) = ---“ j U’(p) exp (p 2 - 2p’p) + ij '2 (p»P )J d P ^ 

where U’(7) is the field solution after reflection from the target, and 
^2 represents the turbulence effect from the target to the receiver. The 
mean intensity at the receiver is then 

<Kp)> = < JU(F) | 2 > = Jj dp 2 dp 2 <U'(p 1 )C , *(p 2 )> 

’ eXp [IT M ~ p 2> - 2p-(p 1 -p 2 ))] w 

• < exp j^ 2 (p,pj) + * 2 *(p,p 2 )] > 
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Through the assumption of a diffuse target, the reflected beam 
suffers a random phase delay from point-to-point over the target, so that 


<u '(p 1 )U , *(p 2 )> = <I( P j)> 6(pj-p 2 ) 


( 6 ) 


Using this in Eq. (5), the mean intensity becomes 

2 + 

<I(p ) > = (YtTl) J dp l < l U(p O | ;> < ex P [iMP.Pl) + V'2 *(P,Pi)J 

where the mean exponential term is unity from considerations of energy 
2 

conservation. The resultant mean intensity at the receiver is then 
simply 


(7) 


<Kp)> 


(-a/ 


dp <|U(p)| 2 > 


( 8 ) 


To complete the solution, we use Eq, (3) with Eq. (8). We note 
that the structure function gives us (r - |ri-r 2 |) 


<exp [iMp.rj) + >h*(p,r 2 )J> = e 

For the focused beam, we then have 


/ r_\ 

" l P o/ 


5/3 


(9) 


|U(p)| 2 > = 


fe ) 2 v/jf 
/ \ 5/3 " 
fc) . 


drjdr 2 exp 


1 + 


2a 2 
o 


ik — — — 
~jT P' (ri-r 2 ) 


( 10 ) 


Carrying out the integration indicated in Eq. (8), involving the Fourier- 
Bessel integral, we have finally 

<^»-£(r u ° 2 f 
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The result for the collini: ted beam Is identical, and in fact could 
be deduced for an arbitrary beam focus (Eq. (1)) through conservation of 
energy: 

<I(p)> = (ife) / dr '''l u (o! 2 > 



0 


Thus the mean intensity at the receiver (illuminator) plane is uniform and 
independent of turbulence level. 

2. Correlation Function of Irradiance 

In order to calculate the correlation function or covariance 
of irradiance, we assume for the present section that the perturbation Green's 
function (wave structure function) is dominated by the phase perturbation 
(phase structure function).'* This will be true for many cases of interest 
and will be relaxed in subsequent sections, where the actual implications 
of the assumption are pointed out. 

The correlation function of the intensity at receiver points pi 
and p 2 is given by 


5 I (Pl,P 2 ) - <Il(Pi)I 2 (p 2 )> = <U(Pi)U*(Pi)U(p~ 2 )U*( P2 )> (13) 


Utilizing the extended Huygens-Fresnel principle, the correlation function 
can be expressed as 


B j(Pl.P2) 



dp I dp 2 dp 3 dp 4 <U(p 1 )U*(p 2 )U(p 3 )U*(p 4 )> 


exp[ik(R 1 - R 2 + R 3 - R 4 )] H( Pl ,p 2 ,p 3 ,p 4 ; Pl ,p 2 ) (14) 


5, R, S. Lawrence and J, W, Strohbehn, "A Survey of Clean-Air Propagation 
Effects Relevant to Optical Communications", Proc, IEEE, 58, 1523, 
October 1970, 
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where H is the fourth order mutual coherence function given by 


H - < exp[«(pj, Pl ) + i|>*(p 2 , Pl ) + ^(p 3 ,p 2 ) + ij/*(p 4 ,p 2 )i > 


and 


R 1 = IPl ~ Pi| 
R 2 = I Pi - P 2 | 
R 3 = |P 2 ~ P3I 
R 4 = I P 2 — P 4 I 


( 15 ) 


Under the assumption of dominant phase perturbations. 


H < exp[ i<. (p j, pi) - is'(P 2 »Pl) + i<J>(p3, P o) - i<^.(p4»P 2 ) ] 


(16) 


After reflection from the diffuse target, the fields are Gaussian 
and spatially incoherent. Therefore, the fields at the target can be 
expressed as 


<U( Pl )U*(p 2 )U(p 3 )U*(n l4 )> 


<u (P1 )U* (P 2 )> <U(p 3 )U*(o 4 )> 


+ <U(p ] )U*(P l4 ) > <U*(p 2 )U(p 3 )> 

= <I(Pi)> <I(p 3 )> 6 ( P ]-p 2 ) 6(^3-^) 

+ <I(Pi)> <I(P 3 )> < 5 ( Pl -^ 4 ) 6(p 3 -p 2 ) 

Utilizing (17) and (14) the correlation function can be expressed as 


4 


B 


+ 


I 



(2^) JJ dp2 dD,+ <I(p 2 )s> < 1 (p 4 ^ > 

dp 2 dp 4 <I(p 2 )> <I(p 4 )> exp[~ (p2 Pl<) (Pl P2 |] 

• H(p 2> P2 » p 4> p 4iPl> P2) 


(17) 


( 18 ) 
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where use has been made of 


R 1' R 2 - 2L p 3 2 “ P 2 2 ~ 2 (Pl-Pj 
R 3 - R M P 3 2 “ P4 3 ~ 2(p 3 - 

The fourth order mutual coherence function (see Appenfix A) in Eq. (18) 

. . , 2,6,7 

is given by 

3 ’(p 2 >P 2 »P4»Pq;P]»P 2 ) = H (P]»P 2 *P 2 »P4;P].P2.P3.P4) 

Pi = P 2 

P 3 = P 4 
P3 " P] 

P 4 = P 2 

2 (D12 “ 33)3 + 3 ) ] 4 + D 2 3 ” D 2 4 + D34) 

Pi = P 2 
P 3 = Pt * 

P 3 = Pi 

P4 = P2 (19) 

where the wave phase structure functions D are given by 


6.. D. L. Fried, "Effects of Atmospheric Turbulence on Static and Tracking 
Optical Heterodyne Receivers", Optical Science Consultants, Technical 
Report TR-027, August 1971. 

7. D. L. Fried, "Atmospheric Modulation Noise in an Optical Heterodyne 
Receiver", IEEE Trans, on Quantum Electronics, QE-3 , 213, June 1967. 
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D ij " p 5/3 


iPj-Pil 873 " I p i-PiI 8/3 | 

------A, 

( PfV ~ <^-^>1 


where 


P„ = (0.545 C 2 Lk 2 ) 
o n 


is the turbuj.jnce-induced coherence 


scale and 


? 2/3 

C n = structure constant of index of refraction (m~ ). 

Using this in (18) and making the change of variables 


P - P2 ~ P 4 , p = Pl - P 2 


2R = p 2 + p 4 

and recognizing that the first term in Q8) equals <I(p I )><I(? 2 )>, the 
covariance for the focused case is given by 


(miQ) i"„r (+)’,■•© 


J o ( C kr >l R ' + l0 J„ ( k r 2 ,R - fl) • 


k r,|R - P-|N ' P 


r ] r 2 drjdr2 ex P “ 


/ r 2 i 2\ , 5 / 3 5/3 

(r l + r 2 / (rj + r 2 ) 


5/3 _ I I p 8 / 3 - p 6 / 3 | 
Ip - p| 


1 | pS/3 _ q8/3 

" 2 - 

Ip + p| 






The covariance for the collimated case is obtained from (21) by replacing 
the dridr 2 integration by 



r 2 drjdr 2 


exp 


(rj 2 + r 2 2 ) 



+ 



5/3 5/3 

r 1 + r 2 

p 5/3 
0 



p 5/3 



P " P 


1 l£ 

2 


8/3 


,8/3 


P + P 



( 22 ) 


In order to further reduce the number of integrations, the zero 
order Bessel functions must be expanded to functionally separate the R and 
p dependence. This can be accomplished by utilizing the following iden¬ 
tities : 


J 

o 




— \ to 

2 _i ) *J 0 (l r ‘ s ) (c r 4) 


cos m <\> 


where 


and 



6=2 

m 


for m ^ 0. 


The 0£ integration is then given by 


(23) 
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/ J o (l r ‘ |R + 2 >) J o (l r2iR ' I 1 ) de R 
0 

00 m /k \ /k \ 

■ 2 ' 2 '- 1 * E „ J „ (l r > R ) 4 (l r * R ) 

ra :i U 

• (r '*!) j. (*- r4) 


Using the Fourier-Bessel integral, 


C2A) 


Jri drjf(r 2 ) J RdRJ ra (^ ri R) 

o o 



and the covariance for the focused case becomes 


= f (r 2 )for r 2 ^0 

2 

_1/L\ f(o) for r 2 =0. 

(25) 



focused 


• exp 




J 


(M t^K] 


exp 

ik — — 2 

- p • P - 

n 5 7 3 _ I |p 8/3 - P 8 /3| 1 | p 8 / 3 _ p 8/3 jj 


p V3 

— 0 

2 ~ " _i 

Ip ~ p l 2 |p + p| 


where lse has been made of 


l 

m =0 


m 

(- 1 ) e 


m 


J 2 
m 


(x) = J q ( 2 x) 


(26) 


(27) 
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and 6 is the Kronecker delta. 
T2 


The variance is then given by 


o 


J 


2 




( 28 ) 


and the normalized variance 



is unity (Eq. 12). 

The same technique yields the covariance for the collimated case 
It is given by 



exp 


ik - - 2 

7- P * P - 


5/3 


p 5/3 _ I |p 8 / 3 - p S/ 3 j 1 ]p 8/ ' 3 - p 8 / 3 

Ip - p| 2 Ip + p| 


and the normalized variance is again unity. 


(29) 
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The above results are based on the assumptions of a diffuse 

target and phase perturbations dominating the turbulence effects. The 

resulting threefold integrations will be numerically evaluated in future 

work. Physical interpretation is not difficult but will be clarified 

in the discussion of further simplifications below. It suffices to point 

out here that the covariance fundamentally involves two scales: the 

turbulence-induced coherence scale (P Q ) and the sperKle scale in the 

absence of turbulence (a and L/ka for the focused and collimated case 

o o 

respectively). 

Alsi, we note that within the present assumptions, the normalized 
variance of intensity is always unity independent of turbulence strength. 
This igrees with a physical model of identical-frequency, randomly 
phaseo oscillators summed to -epresent any given point in the receiver 
field: the model applies regardless of whether target speckle (o q or 

L/ka Q ) or "atmospheric speckle" (p Q ) dominates. 

3. Mutual Coherence Function 

The mutual coherence function (MCF) may be very important in 
analyzing the operation of a coherent optical adaptive system, and can 
be readily derived given the assumption of a diffuse target but without 
assuming dominance of the phase perturbation term. We write 


2 

r(pi,P2) = mu dpjdp 2 < U(pi)U*(p 2 ) > exp jikjRi(Pj,p^) 

- R 2 (p 2 ,p 2 )]|< exp ^ 2 (Pl*Pl) + 4'2 *(P2>P2 )] > 

where Ri(Pj »Pi), R 2 (P2>P2) are'the distances from pj to pj and p 2 
to p 2 respectively. 


(30) 


By the Fresnel approximation 


R 1(Pj >Pi) 


R 2(P2 > P2) 




P] ‘ Pi ~ P 2 ‘ P2 

L 


( 31 ) 


- 15 - 



Finally, from (30) and (31), 


r( ^>' (i-ir) 


exp 


i k (pf-ps T| 
2L 


fC - 

jj d Pl d P2 <U(f-i)U”(P 2 ) > 




< exp 


X P . 2 CPI, P i ) + (p 2 > P? ) 


(32) 


Since the wave is incoherent after reflection from the diffuse 
target, the coherence function at that plane can again be represented by 
the Dirac delta function as given in ( 6 ). Using this in (32), r(p lt p 2 ) can 
be simplified to 



ik (P?-p|)' 


2L 


I r _ 

j de <I(P)> exp At (p,-p 2 ).-| 


(3’) 


In the absence of turbulence, this equation is entirely identical to the 
Van Cittert-Zernike theorem of coherence theory , 9 which is identical to a 

result obtained by Goodman for the mutual coherence function of a pulsed 
optical radar . 10 

To complete the solution, we utilize the mean intensity at the 
target. For the focused case: 



We thus have 


9. M. Born and E. Wolf, P rinciples of Optics . Perg&mon Press, New York, 
1975. 

10. J. W. Goodman, "Some Effects of Target-Induced Scintillation on 
Optical Radar Performance", Proc. IEEE, 53, 1688, 1965. 
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r <pi*P 2 > ■ (ik) (l) Kl i-J dp / r dr J o (l 


Ad 

• e o 


r 

2 \P 



5/3 


exp 


ik- 

r p,p ' 



5/3*1 


exp 


ik(pT-pi) 

2 L 


/ \ 4 2 a 2 

Ul) ' U ol -f- 


CD CO 


/ 


pd P J* rdr J Q (£ P.) J D (£ PP. 


* e 


4a 2 ' o 

o 



5/3 , ,5/3 



ik 2 

2 L 


(Pl"P*) 


(34) 


From the Fourier-Bessel integral formula. 


CO 

/ 


PJ 0 (l rp ) J o (c dp ■ (I) ^ 
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(35) 


Equation (34) can then be simplified and it becomes 


r(pi.P2> - h (r^ l«, l 2 -S- e 


2 (£- 
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ik 
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5/3 
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= <I(p)> e M < 36 > 

Using <I(p)> for the collimated case in (33) and simplifying yields 
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- p 2 


fc)‘*(Sf - fef 


ik 


( 37 ) 


2L ( p ?-°j) 

r(Pi.P2) = <I (p)> e 1 

Collimated 

These results for the MCF will be used further below. 

It may be nct*d that the MCF's of Eq. (36) and (37) imply a "white” 
or constant spatial power spectrum for the (complex) amplitude. This is 
of course an idealization resulting from assuming delta-function rather 
than wavelength-sized phase correlation for the field upon reflection off 
the target. The more interesting spectrum, however, is that of the irra- 
diance, as discussed below. 

4. Probability Distribution 

In order to formulate the probability distribution for the scin¬ 
tillating energy at the receiver, we evaluate the n th moment of the intensity 
in terms of <I(p)>. We again assume that the phase perturbations are the 
dominant turbulence effect, and write the second moment as 


<i2/~n> _ / j<_\ 4 fCff - _ _ _ 

?) \ 2l,L /J ijj dp l dp 2 dp 3dP 4 

• <U(pi)U*(P 2 )U(p 3 )U*(p 4 )> exp - p | 

+ - pg - 2(p~j~P2) • P^^-p,) j 


<ex P i(^(Pi,p) - <Kp 


2 .p) + * ;p 3 ,p) - <kp 4 ,p)) 


(39) 


Since the fieldd after reflection from the diffuse target are jointly 
Gaussian and uncorrelated, 
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<U(p 1 )U*(p 2 )U(p 3 )U*(p lt )> = <U( Pi )U*(p 2 )> <U(p 3 )U*(p 4 )> 

+ <U(p ] )U*(p 1< )> <-U*(p 2 )U(p 3 )> 

I(P])6(p 1 -p 2 ) I (p 3 )6(p 3 -p 4) + I(pj)6(pi-p 4 )l(p 3 )6(p 3 -p 2 ) 

Using (40) in (38) yields 


3~P2-> (40) 


-I 2 (p)> = 2 


a // 


dp 2 rfp 4 < I(P 2 )> <I(p 4 )> 


(41) 


Since 


<P)> ' (fk) / dp 2 <I<P 2 )> 


<I<P)> ’ (■&) 
the second moment becomes 


< I (p) = 2< I (p) >2 


(42) 


Similarly, 1c may be shown that the n moment is given by 

<l"(p)> - N! <I(p)> N (43) 

The probability density function for the intensity therefore is exponential 
and 

a 

<I> 

(44) 

where a is greater than or equal to zero. It is thus concluded that the 
field and amplitude at the receiver are normally and Rayleigh distributed 
respectively, given the assumption that phase perturbations dominate the 
turbulence effects. 


v Pj (a) = 
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5. 


Simplification for Weak and Strong Turbulence 

Since, as shown above, the fields at the receiver are gaussian 
and spatially "white", it is tempting to assume that the receiver fields 
are also jointly gaussian. This turns out to be a good approximation in 
many situations, and in this section the implications and conditions for 
validity of this added assumption are explored. This leads to a simple, 
straightforward interpretation of the terms in the covariance of intensity. 

The jointly Gaussian assumption yields 


B 


I 


<U(p))lT ;: (pj )> < U(p 2 )b A (p 2 ) > + <u (Pl) u, ’‘ (P2) ><U * (P’ (P2) > 


= <l(p])> < I (P2) > + !^(p1,P 2 )I 
It follows that the covariance of intensity is given by 


(45) 


C T (p 1 > P2 ) = b I (pi,p 2 ) - <1 (p])><-I (p 2 )> = I f (p ]»P 2 ) I (46) 

Finally, utilizing the mutual coherence result (Eq. 36) the 
normalized covariance function of irradiance for the focused case can thus 
be written: 



(47) 


focused 

where the normalized variance is unity as before. 

For the collimated case, the same variance is obtained, and the normalized 
covariance is 


C T (p) 
N 



collimated 


(48) 
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The covariance scale lengths are obvious from these results. 

Either the "atmospheric speckle" (p Q ) or the "target speckle" (speckle 
in the absence of turbulence) will predominate, depending upon which is 
smaller (strong and weak turbulence respectively). We point out in passing 
that a third covariance scale (L/k) may also enter, but this scale is lost 
within the present assumption of dominant phase perturbations (see Sec¬ 
tions II-B and II-E below). 

We note that the spatial power spectrum of inadiance may be 
readily obtained by transforming Eqs.(47,48),, However, a more important 
quantity in the operation of e.g. an adaptive optical system may be the 
temporal spectrum. This spectrum, which will be derived in a later sec¬ 
tion, depends only on the atmospheric speckle term; the target speckle 
field will not translate with the transverse wind. 

We now explore the conditions for validity of the jointly gaussian 
assumption. The simple multiplicative terms for the covariance scales in 
Eqs. (47, 48) are replaced bv a more complicated interrelationship in 
Eqs. (26,29). It may therefore be surmised that the jointly gaussian 
assumption is valid under conditions of weak and strong turbulence, when 
target and atmospheric speckle terms respectively predominate, but that 
the jointly gaussian description is not correct in the range of inter¬ 
mediate turbulence effects when both scales are important and interact. 

We now show that this is indeed the case, 

Weak Turbulence 

For the weak turbulence case, p >>/L/k and the term 



(49) 

in (26) and (29), The covariances then become identical to those derived 
using the jointly gaussian assumption. 
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Strong Turbulence 

In the strong turbulence case, p o «/L/k, which corresponds to 
multiple scattering, or "saturation of scintillations" for a point source. 
Let us consider Eq. (26) with p^-K). The only interesting range of the 
argument (p) is 0 <p<p . The Bessel term (^pJ in polar coordinates) is 
appreciable only for p#G, i.e. P>>P Q or p>>p; and because of the latter 
condition the final bracket in the equation is zero. Hence the condition 
(49) is again obtained, and the covariances again become identical to 
those derived using the jointly gaussian assumption. The atmosphere has 
"decoherentized" the field in a manner similar to that of the diffuse 
reflector itself, 

6. Time-Lagged Covariance and Temporal Spectrum of Scintillations 

For this development we assume that the fields are Jointly gaussian 
at the receiver and consequently 

Vpi.pz.t' - |r ( p 1 ,p 2>T) | 2 ( 50, 


Using the extended Huygens-Fresnel principle 


ik(p1 2 -p 2 2 ) 


/ — — \ 


" 2L 

cc 

I \Pl»P2>V 

‘ ( r 'h 

) 6 

Jj dpjdp 2 

• <u(pi,o) 

'K* (p 2 , 

t) > CX p (j£ 

(Pi 2 " P2 2 ~ 2 Pl ‘ Pi 

- 2p 2 P 2 ) 

|J <ex P 

V'2(pl>Pl»°) 

+ *2*(P2>P2> T ) > 


(51) 


Due to the diffuse target 

<U (P1,0)U* v'p 2 ,t)> = (p"j , 0)U*(p" 2 , t) > 5 (p*! -p" 2 ) 


(52) 




Using (52) in (51) and utilizing the extended Huygens-Fresnel principle 
to express the fields incident on the target, the time delayed mutual 
coherence function at the receiver becomes for the focused case 



(53) 


The first ensemble average in (53) corresponds to the time delayed 
mutual coherence function for spherical waves originating at two points 
rj and r 2 in the transmitter plane and propagating to a single point pj 
in the target plane. The second ensemble average corresponds to the 
time delayed mutual coherence function for a spherical wave originating 
at the point pj in the target plane and propagating to two points pj and 
P 2 in the receiver plane. Performing the integrations in (53) it becomes 


/ . 2 / . 2 2 00 
r( p-<>' (si) (t) V -f-/* / 


rdr J o(r pr ) 


exp 


r* .k- 

-- - i- p-p 


4a 


F(r,0,t)F(0,p,T) 
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-£G) I’.i’V- 


4a 11 
o 


+ 51 (Pi 2 -P2 2 ) 


F (P»0,r)F(0,p,x) (54) 


where the time delayed mutual coherence functions F(p,0,t) and F(0,p,t) 




can be obtained from the mutual coherence function for a spherical wave 
by invoking Taylor's Hypothesis. 5 

This mutual coherence function is given by 2 


^O.P) - exp -2.91 Lk : 


/ C n 


('■’) l w P + (l-w)rl 1 dw 


(55) 


where r is the transmitter aperture vector r 1 -r 2 ,p is the target aperture 
vector Pi-p 2 and w is the distance from the source to the field point 
normalized by the total path length L. For the uniform turbulence case, 
we note that k 2 LC^ 2 %p o ~ 5 / 3 . The time delayed mutual coherence function 11 
can be obtained from (55) by replacing wpT by wp - V(w)t, where V(w) is 
the transverse wind: 


F (r,p,T) = ex 


<P 


-2.91 Lk 2 


J c n 2 (w)|wp - V(w )t + n-v;)r| 5/3 dw 


and consequently 


(56) 


F(p,0,t) = ex 


exp 


ro target 


and 


F(0,p,t) = exp 
From target 


1 

2.91'Lk 2 f r 2 

0 

] 

/* 


_ _ 5/3 

C n z (w)|(l-w)p - V(w)x| dw 


•] 


' (57) 


2.91 U 2 j C n 2 (z)| zp - V(z)t| 5/3 d 


(58) 


where z = (1-w). 

The normalized, time delayed covariance function for the focused 
case and an arbitrary distribution of C n 2 along the path is thus given by 

,2 


_ _p: 


1 

r 


C L . (p-,T)--e 2« 0 2 „p[- 3 . 82 j c n *(»,| ( i^p.v(. )T | 5/3 


N 

focused 


dw 


11. R.S. Lawrence, G.R.Gchs, and S. F. Clifford, "Use of Scintillations^ 
Measure Average Wind Across a Light Rean". Applied Optics,11, 239 
February 1972. —* 
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For the collimated case, a further multiplicative factor 



completes the expression. 


The temporal socctral density of irradiance can be obtained from 
(59). Letting p = 0 and assuming uniform turbulence and crosswind, (59) 
becomes 


C T (0, t ) 
N 


-5.82 Lfc 2 C 2 |v< 5/3 t 5/3 

n ! i T 


-10.67 

e 



(60) 


Taking the Fourier transform of (60) yields the spectral density: 


Sj (to) 


= 2 e 


I 


-10.67 


5/3 


5/3 


Cos( u t) d t 


5 |V| (10.67) 3/5 


rn/ rx X 2 r(9/5) 
r(3 /' ,) - -T/5 + 


2 ! 


(10.67) 




2(n-1)! o (n-1)/5 

V } (10.67) 


(61) 


where 



M 


The normalized spectral density is plotted /ersus the parameter (x) in 
Figure 2. 
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Figure 2. flonnalized temporal spectrum of irradiance scintillations for 

a focused beam, diffuse target, and jointly-gaussian field. The 
frequency (w) and spectrum (Sj) are normalized by the uniform 
transverse wind speed |v"|and the turbulence coherence radius p as shown 



II.B.The General Covariance Function 


In the previous section, the mean irradiance and mutual coherence 
function were derived with the simple assumption of a diffuse target. How¬ 
ever, the analysis of the intensity covariance was based on the added 
assumption of neglecting the amplitude term in the perturbation Green's 
function. This latter assumption is not always valid, as can be seen by 
considering the case where the illuminated target spot in small: the scin¬ 
tillations at the receiver should then approximate those of point source, 
with log normal statistics and a covariance scale on the order (L/k)* 5 . 

In this section we include the amplitude perturbation term and examine 
the effect on the probability distribution of the irradiance. We then 
derive the general covariance function, with physical interpretation. The 
actual conditions for validity of the simpler result of Sec. II.A will also 
be discussed, with further elaboration in Sec. II.E. 

1. Moments and Probability Distribution of Irradiance 

The second moment of irradiance can be written from Eqs. (14) and 
(17) and the mutual coherence function (15). We now generalize the latter 
to include amplitude perturbation terms (Appendix A), resulting in 

<I 2 (P)> = 2 (-“r) JJ dp 2 dF 4 <I(p 2 )><I(p 4 )> e 4Cx ( l p 2~ p 4l) (61) 

where is the point-source log amplitude covariance function given in the 
first-order theory ’ * ^ by 

“ L 

C X <|P 2 -?,|) - j u*(u>J o (l£2^d us ) Mn»[2?f fipil ] d , du (62) 

0 0 

and $ is the Kolmogorov spectrum^ 


®(u) = 0.033 C z u 
n 


2 "ll/- 


(63) 


!!_!!!!!! Jl r the present that V* whlch wil1 be tru * 5 for both wea k 


12 ’ Mon T N»M 8 kl i’ ~ ,e ~ E f f e ? tS ° f the Turbule nt Atmosphere on Wave Proofs 
,, ~—J' National Technical Information Service (0TT-68-5O464), 1971, P * 

* ee an< * H arp, "Weak Scattering in Random Media with AnnH 

cations to Remote Probing", Proc. IEEE, 57, 375, April 1969. PP 
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and saturated scintillations (weak and strong turbulence scattering respec- 

tively). It will not be true for the intermediate case, however, and we 

will lift the restriction in calculating the covariance below. Also, we 

point out that the function C has been derived phenomenologically for the 

X .I* I/: 

saturated or strong scattering case by Yura and Clifford so that in 

principle we are not limited to the first-order expression (62), 

We generalize (61) to the nth moment: 


<!"(?)> - „* 


2 n 



<I(P ± )> 



2 


n 


E 

i*J 



We write this as 


(64) 


<l n (p)> = n!<I(p)> n F 

n 

where 

F], “ 1 

F 2 JJdp 2 dp4<I(p 2 )><I(^ 4 )>e 4C X ( l p2 ’' p ^ ) 

>p 2 <I(p 2 )>] 2 
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• ___ _ s 

dpidpy , . .dp" 

n 

[j, <I(o i )> ] 

exp 

[2 j c/i^-pjD] 


[ Jdp<I(p)> ] 

n 


Our assumption on C^ yields 
4 C X (IP 2 -P 4 I) =1+4C (|P 2 -P4 I ) 

^ A. 

We thus simplify F 2 to 


(65) 


( 66 ) 


(67) 


14. 


15. 


16 . 


4ffdp 2 dp4<I(p 2 )><I(p 4 )>C (|p 2 -p 4 |) 
Fo - 1+ — -X_ 


[$dp 2 <I(p 2 )>] 


( 68 ) 


H. T. Yura, "Physical Model for Strong Optical^-Amplitude Fluctuations 
in a Turbulent Medium", JOSA, 64, 59, January 1974, 

S. F. Clifford, G. R. Ochs, and R, S, Lawrence, "Saturation of Optical 
Scintillation by Strong Turbulence", JOSA, 64, 148, February !'74, 

S. F. Clifford and H. T, Yura, "Equivalence of Two Theories of »trong 
Optical Scintillation", JOSA, 64, 1641, December 1974. 
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For the denominator for the focused case, we write from Eq. (10): 


:i<p)> '(l) l“J 2 -f-J rdrJ 0 (? r ) 


r 2 / r v*/ 3 

m 2 "UJ 


which integrates to 


<I(p)>dp = 2ir|U I 2 

o 2 


The integral in the numerator of (68) is straightforward but laborious, 
We state here the result for F 2 : 


u 2 s 2 /us \ 5 / 3 

L " 2a 2 k 2 - \*r) 

F 2 - l+(4u) 2 k 2 J uduj" dse ° $(u)sin 2 


0 0 


where the approximation relates to the assumption (67) on Cx. 
We then have the second moment from Eq. (65): 

<I 2 (p)> = 2<I(p)> 2 F 2 

and the n^ moment is 


:i“> - n!<I> n jl+ (F 2 -l) j 


Hence nonunity F 2 represents the departure from an exponential intensity 
distribution. 

We finally derive the probability distribution using 


< e > * l-s<I>+ 


'I 2 > s 3 <I 3 : 


+ ... 


We let <I> = I and find the probability 


VD-i-s 0 i+ <,,_,)[ i-|i + £,] 

u o o J 

so that we have a first-order correction to the exponential distribution. 
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We note that Fo (and F ) will be unity for weak turbulence and for 

n 

very strong turbulence (p^-K)) , Also, F 2 will be unity for a small focused 

source (a ), which physically relates to a large target spot. Conversely, 

0 

F 2 will depart from unity and therefore show the effect of the amplitude 
perturbation term for the case of intermediate turbulence and a large 
source (small target spot). Physically this is the case of a quasi-point- 
source attempting to scintillate in the usual manner for a point source 
in the first-order theory, 5 but nevertheless interacting with the speckles 
created by the diffuse target. 

The exponential distribution applies for the weak turbulence case 
simply because the behavior is that of a diffuse source in a vacuum 
(speckles). It applies to a large target-spot, implying that the phase- 
perturbation-dominance assumption of Section II.A is then applicable, i.e. 
that the target (vacuum) speckle mechanism dominates. It applies for 
strong turbulence because the field again has the nature of that from a 
diffuse source (atmospheric speckle p^). These considerations will be 
clarified below and in Section II.E. 

2. Covariance and Variance of Intensity 

We now derive the general covariance function, dropping the requirement 
that C <<1. 

x ,, ON 
The covariance is the second tenr of Eq. (18): 

/ u \4ff „ _ _ T^(Pl-P2)(Pl“P2) __ 

c I (Pl»P;>) “ \2 ^l) JJ dp 2 dp 4 <I(p2) > ' I ( p 4) > e H(p ,p) (76) 


where 


P 51 P2rP4 and P = Pl"P2 • 

The full coherence function replacing Eq. (19) is (Appendix A) 


H(p,p) =* 



[ p 8/3_ p 8/3| 

p 5 /3 
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(77) 


2C (p,p)+ 
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2C (p,-p) I 

A ' 
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where 


17 

1 

<yp,p) = 0.132^k2C n 2 L j dtj du u' 8/3 sin 2 

o o 

•J o [u|tp + (l-t)pl] (78) 

The latter expression (78) assumes that multiple scattering (saturation) 
does not apply but can be modified for such a case,^ The irradiance 
at the target <I(^)> is given by (69) for the focused case and with an 
additional (ka Q /2L) 2 term in the exponent for the collimated case as before. 

We again omit tedious algebra and integrations and state the result 
for the focused case: 



J r 2 dr 2 l 

o 





ik- 

-r p *p_ 

H(p,p) 


(79) 


The variance is given by 


Cj(o) = 


- (h) (tf l»J* 


2 a 

(/-) J 


r 2 dr 2 


2a 2 
o 


2r 2 5/3 

p 5 / i 

0 


• J »d 0 J o (ir 2 p) e‘* C X <o) 


(80) 


These^general results involve five-fold integrals (more in the saturated 

case ) and have not been numerically evaluated. A comparison with Eqs. 

(26) and (28) show that the amplitude perturbation term simply introduces 

the _^ CermS ln Eq * ^ 77 ^’ i,e * the lo 8 amplitude covariance for a point 

17. T.Wang,S.F.Clifford, and G.R.Ochs, "Wind and Refractive Turbulence Sens¬ 
ing Using Crossed Laser Beams", Applied Optics, L3,2602, November 1974. 
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source. A qualitative interpretation will he given below. 

In order to compare these results with those of Section II.B.l, we 
again let C^<<1, and we find that (80) and (63) yield 


.1 


» 1 + 


0.528ir 2 k i: C^L ^ dtj* du u exp J~2 L(l-t)J‘ 


- ^5/3 [f L(l-t>] 5/3 sin* [ u2 ^^' :t) ] (81) 

0 

This is consistent with Eq. (72), and again shows the departure from the 

unity normalized variance obtained when C cannot be taken as zero. 

X 

3. Further Remarks on Probability Distributions and Limiting Conditions 

A systematic, quantitative description of parameter realms will be 
given in Section II.E. However, we can indicate here the types of considera¬ 
tions involved. Let us write the general coherence function (77) as 


- 9 ^- 2 —\ 


H(p,p) = e 


H'(p.p) 


where 


H'(p,p) = e 


-e ) 5/3 ♦ 


,8/3_ p e/3| 

P o 3/3 


f^r + wr> 2c * (p ’ p)+ 2 V p> - p) 


From the preceding discussions, we know that non-negligible terms relate 
to a nongaussian field (irradiance is not exponentially distributed, and 
that H ,,= l corresponds to a jointly gaussian field. The amplitude perturba¬ 
tion effects (C^) will be influential under the following conditions (see 
Eqs. (79) and (83)): 

A. Turbulence is moderate or weak (i.e. p >(L/k) , and 

1* 0 

B. The source is large (a Q >(L/k) ), 

These conditions can be reasoned from Eq. (79) by considering the large 
range of t 2 in the exponential, leading to a small range of p from the 
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Bessell function, and the .o,-responding influence of C in (83) on the 

X 

result of the p integration. Physically, this corresponds to nonsaturated 
turbulence, with a near-field target and small spot acting as a quasi¬ 
point-source . 

Conversely, if o q is small, the range of T 2 is limited, and that of 

p expanded such that does not make an important contribution to the 

p integral. This means that a large target spot results in dominance of 

the phase perturbation terms. In fact, if p »u , H' can be taken as 

o o 

unity, the field is jointly gaussian, and the covariance is given by Eq. (47). 
Finally, if turbulence is strong ( Pq <(L/ k)**), which corresponds to the 
multiple scattering realm, then with a large or small source the p-integral 
will be controlled by the phase terms in (83). Thus the phase perturbation 
dominates in strong scattering; the field will be jointly gaussian except 
in the transition region P q Jfc c^, where it is simply gaussian. 

4. Further Discussion of the Covariance Function and Scales 

The qualitative behavior of the covariance curve, with the attendant 
covariance scales or scintillation patch size, can be described as follows. 
Suppose that we have a large, focused transmitter (near-field target), 
with weak turbulence ( a 0 ,P 0 > (L/k)* 5 ). The scintillations will be essentially 
those of a point source (unsaturated): 5 


0.124 C 2 k 7 / 6 L u / 6 


n 


(84) 


for the log amplitude variance. The covariance curve for irradiance will 
be as^shown in Figure 3a, where the covariance scale is proportional to 
(L/k) . Now let the turbulence strength increase (p decrease) until p 
is comparable to (L/k) 1 *; the log amplitude variance peaks at approximately 
0.6 before saturating, 5 and a new covariance scale (p Q ) appears as concep¬ 
tualized in Figure 3b. Finally, let the turbulence become stronger yet; 
the normalized irradiance variance is unity and the covariance curve is 
as given in Figure 3c. The third possible covariance scale, a itself 

7 q * 

will only predominate when o q is small (far-field target). In the latter 
case only, the speckles will not be transited by the transverse wind. 
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Figure 3. Conceptual normalized function of irradiance for a large, focused 
source and varying turbulence strength, o 2 denotes the log 
amplitude variance for a point source. * 

A. Weak turbulence (p 0 >>(L/k)^) 

B. Transition region ip 0 %(L/k)^ ^ 

C. Strong turbulence or saturation (p <<(L/k) . 



This will also be true for a large target spot arising from defocusing or 
collimation, such that the associated speckles are small, 

II.C.Target Glintf 

As a preliminary examination of the effects of target structure, we 
present in this section an analysis of the mean receiver irradiance in the 
presence of "glints" (specular reflectors), We assume that the target is 
otherwise diffuse, and that phase perturbations predominate. The transmit¬ 
ter is again a TEM^ focused or collimated source, 

We again start with the field at the receiver as described by Eqs. (1,2). 
The field reflected at the target is represented as 


"<0> ■ "diffuse ( 0> + "glint "g^ 


where 
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U= Z U . . (p ) a exp 

glint , incident M m m 

b m=l 


p - position of m th glint 

® til 

Ap = width of m glint 

^ th 

a = complex strength of m glint 
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-(p-p ) 2 
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Ap. 
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Then the irradiance is 


(85) 


«P> ' (2^L) 2 Xfl"d I ^ + U g<P1 |> d *<P2>+U‘<P2>] 

’ exp { 2L [^ p l - P^ 2 'fP 2 - P) 2 ] + *2 (Pi ,P)+1<* (P 2 .P)) 1 d Pj - dp 2 

* [";<P _ 2) + "*(P2)] 

. exp [p?-pf -2(p]-p 2 )p] + li 2 ( p 2 »P)'H , 2 (p 2 ,p) ^ dpj ,dp 2 (86) 


We then simplify the field at the target: 


-33- 


[ u d ( p i ) +u g ( p :)] [u d (P2)+U g (p 2 )] > 

= <U d (p 1 )U*(^ 2 )> 4 <U d (p 1 )U*(p 2 )> + <U g (pj)U*(p 2 )> 

+ <U e (Pl)U(P2) > 

O & 


Since 


<U d (pl)U e (p 2 )> = < U d (p I )><U*(p 2 )> = 0 


— * _ 


<u (p l )u d (p 2 )> = <U (Pi)><U (p 2 )> = 0 


<U d( p i)U*( p 2 ) > = <I(p J )> 6 (pj-p 2 ) 


^2 Cpi ,p)+i^ 2 (P2,p) 


I /3 


we have 


<I(p)> = <I(p)> d 4- <I(p)> 

8 

where <I(p)> d i s the diffuse term analyzed in Section II,A 
the glint term: 


and <1(p)> is 
g 


<I(P) Y ( 2 ^) W dp l d P2 <U g (Pl)U*(p 2 )> 

• exp {f [pi-p^Cpj-p.xj]- i^ 3 5/3 } 

We now describe the complex correlation of the glint field 
target: 


at the 
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m 2 


rdr J 


where for the focused case 

-(r) 2 l“J 2 r } 

o 

The corresponding expression for the collimated case has an additional 
(ka Q /2L) 2 term in the exponent. The coherence function can be readily 


_r£ (r \ 5/3 
/k \ ' *ct 2 ' 

(rv) e 0 0 
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shown to be 


<U(p )U*(p )> - i- iu I 2 ^2 

m l 2 tt \ L / 1 o 1 2 
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4a 2 2L 1,11 m 2 
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where the mutual coherence function F is 


F ( r >P -p m ) 
m^ m 2 


= F(r,p) = e 


2.91 C 2 L k 2 T 3 p8/3 , r 8/3| I 

0 n -jr --I- 

1 8 p-r[ J 
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Finally, the mean receiver irradiance (diffraction pattern) due to glints 
is expressed from (90,91) as 




/ k \ 2 M M 

p)>8 ' fe ) 4-, 4-1 V^V'V* 

JJdPld 02 exp | IK [4-P2-2(?,-p 2 >p] /p\ 

P ">1 Ap m ? V p o/ 


(95) 


where <u (P mi ) u *(P m2 )> is given by Eqs. (93). We now show specific 

1. Single glint (M = 1) 

From Eq. (95) we have 


cases. 


I(P)> 8 ‘( 2 s ;) 2 <I(p n, )> l«J 2 


2 2 

ff f ik r 2 2 - - -1 Pl+P 2 -2p .(p‘ 1 +p‘ 2 )+2p 2 5/3) 

JJ d Pl d P2 exp [_ p l-P 2 - 2 (Pi-p 2 ).pJ-_ _ ' 

P m vo' i 


With the following change of variables: 
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we have 
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The first integration yields 


(98) 
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which with (98) gives 


(99) 


<I(p)> s “(l) 


<I(p ) > |a | 
m m 
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(r) 2 Hy 


5/3 

dp 


( 100 ) 


This may be usefully normalized by setting p - p . 

m 

The result (100) may be interpreted as follows. The maximum occurs 

at p = p m> as expected. The scale of the diffraction pattern (|p - p" |) 

is seen from the Bessel term to be L/kp , where p is determined from 

max max 

the exponential as the smallest of Ap , L/kAp , and p respectively. The 

m m o 

corresponding I(p) scales are L/kAp (diffraction from the glint), Ap 

m m 

(geometric reflection term), and p (effective coherent glint size in 

o 

strong turbulence); the largest will predominate in each case. 

2. Arbitrary Number of Glints (M) 


This involves lengthy but straightforward calculation. The result 


is 


M . M 

<I(p)> = E <I(p)> + 4 E F (I) 

g i=1 2 iflj i,j 


( 101 ) 


where 
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Each pair of glints results in 
irra'Hance locate at 
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a cross-terra bright spot in the receiver 
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This is illustrated in Figure 4. 

II.D.Incoherent Illumination 

"°" depart fr0m Che C8Se 0f TEM 00 laser illumination and briefly 

multi “i / n0n "'’ 0n ° Chr0n ' 8,:ic case - » h id'> -ight include that of a highly 
multimode laser, in principle, all of the expressions of the pre. 

ceding sections can be generated „ith the utilization of a 
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Fourier transform or partial-coherence representation of the source. As 
the source becomes monochromatic, the first-order effect will be to wash 
out the speckles (a Q and p Q terms in Eq. 47). This will obviate the influence 
of the phase-perturbation terms (Secs. II.A,B) and the remaining scintil¬ 
lations will be due to the amplitude-perturbation term and to geometrical 
considerations, as will be seen below. As the spectral width of the source 
grows further, scintillations e.g. for a point source will be further 
smoothed by "spectral diversity", but this is a weak effect requiring a 
substantial fractional spectral width. 

The situation may be simply analyzed if we assume a completely in¬ 
coherent but quasi-monochromatic target spot. We start with the intensity 
correlation function at the receiver from Eq. (14), and utilize delta- 
function correlation of U(p^) with U*(p 2 ) and also U(p 3 ) with U*(p 4 ) to 
write immediately 


<i(pi)Kp2) > ** (i^l) j*J* dPldp3 l U (Pl')I 2 I u (p3 


exp | 4 Cx(p 3 ~Pi,p)| 


(105) 


As expected, tne target-speckle (a Q ) and atmospheric speckle (p Q ) terms 
disappear, leaving the log amplitude covariance function. Since this 
function is unevaluated within the Huygens-Fresnel approach, recourse is 
again necessary to first-order or saturation theories as discussed in 
Section II.B. It is interesting to note however that Eq. (105) is con¬ 
sistent with Eq. (8) of Ref. 19, which was obtained through the more 

. 13 

complicated phase-screen approach. 

In the case of a small target spot (diameter<(L/k) ) and unsaturated 
scintillations, Eq. (105) is a statement of point-source behavior. For a 
larger target spot, it appears that the covariance scale will be smeared 
by geometrical effects, and this can be directly deduced from the spatial- 
spectrum representation of Ref. 19, Eq. (9). (Note however that the quan- 

18. D.L.Fried, "Spectral and Angular Covariance of Scintillation for Propa¬ 
gation in a Randomly Inhomogeneous Medium".Applied Optics,1_0,721,Apr.1971 

19. S.F. Clifford,G.R.Ochs,and Ting-i Wang,"Optical Wind Sensing by Observing 
the Scintillations of a Random Scene".Applied Optics ,14,2844, Dec . 1975. 


-39- 



tity q in that reference is necessarily real, contrary to the notation 
there. Also, q and y in that reference may be identified with |u| or 
I so that the calculated covariance is closely related to but not 
identical to that for I.) It is known from other analyses 20 that the 
scintillation variance is reduced or "smoothed" by a featureless, extended 
incoherent source, specifically 


^ - [D/(L/k)] - ?/3 (K 

where D is the source diameter. This is also the "aperture smoothing 
factor" for an incoherent receiver in unsaturated scintillations. 20 

In the case of saturated scintillations ( Pq < (L/k)S , the covariance 
function in Eq. 105 must be suitably modified. The behavior of the 
irradiance correlation function, in particular its scale size, is not well 

understood for this case of incoherent illumination and strong turbulence 
scattering. 


II.E.Comparison of Alternative Approaches; 
Behavior vs. Parameter Realms 


Variance and Covariance 


The basic method of analysis of this report is the extended Huygens- 
Fresnel principle, which is uniquely suited to expressing problems of 
arbitrary sources in terms of point-source quantities. It also has the 
virtue of inherently retaining the "saturation" or multiple scat'- -ring 
effect including the explicit influence of the quantity Pq e.g tor 
an extended, coherent source. 

The most elegant alternative method is the phase-screen approach, 13 » 19 
a first-order (single-scattering) analysis which also handles an arbitrary 
source (and receiver) geometry but with severe limitations on the strength 
of scattering or phas- distortion for extended coherent sources. It appears 
that this approach can be made to handle saturation effects properly 
through modification of a log amplitude covariance function, 14 " 17 although 

general applicability e.g. to p^effects over a coherent source has not yet 
been clearly demonstrated. 

20. R. F. Lutomirski, et al, Degradation of Laser Systems by Atmospheric 
Turbulence, Rand Report R-1171-ARPA/RC, June 1973. 
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The two approaches are basically complementary when properly com¬ 
pared. The Huygens-Fresnel approach is much simpler in its handling of 
extended monochromatic sources (including the diffuse target case), and 
also simpler in its setting up of the extended, incoherent case. However, 
in general it leaves unevaluated the point—source log—amplitude covariance, 
which must then be obtained from the first-order theory and saturation 
arguments described above. Also, source structure may be advantageously 
represented in terms of its spatial spectrum, as shown in Ref. 19. 

A very useful "hybrid" approach would be to utilize the saturation 
C x function in the Huygens-Fresnel results for various sources. Also, 
it is of interest to explore the behavior of that function for nonmonochro- 
matic sources: how sensitive is the influence of to spectral width, 
i.e. is "washed out" as readily as interferometric phenomena such as 
speckles? 

We are now in a position to summarize the variance and covariance 
behavior, for a variety of sources, in a particularly revealing manner. 

We will identify the fundamental independent variables and all possible 
permutations of their mutual inequalities (asymptotic parameter realms), 
and attempt to show the limiting value of the variances and covariance 
scales for each such realm, and for each source of interest. In order to 
do this we will utilize results from phase-screen or first-order theory 
and from Huygens-Fresnel analyses when respectively appropriate. This 
type of summary has the advantage of clear physical interpretation and 
lacks only the quantitative details of the transition between parameter 
realms. Also, in some cases it permits the deduction of e.g. covariance 
scales through dimensional analysis or the requirement of consistency with 
adjacent parameter realms. 

We present these results in tabular form, with associated remarks 
below. In all cases, the right hand column represents the saturated or 
multiple scattering realm as opposed to the single-scattering realm in 
the left-hand column. Also, "near-field" and "far-field" refer to the 
position of the (point) on-axis receiver relative to the source size. The 
equalities shown omit numerical constants. 


- 41 - 






The symbols are defined as follows: 

Path length: 

Wavenumber: 

Coherence radius: 

Gaussian-beam transmitter radius: 

Incoherent transmitter diameter: 

Covariance scale (scintillation patch 
size) at receiver: 

Beam size at receiver: 

Log amplitude variance: 

Normalized irradiance variance: 


L 

k 

p o 

a 

o 


b 


P 


a 



Note: for log normal scintillations a 2 = e * -1, 

Table I: X N 

Focused TEMqq beam. Realm (a) is the first-order point-source realm, 

while (b) is the "transmitter aperture smoothing" realm within the first- 
, , 24-27 _ J 

order theory. It is interesting to note that the latter has the same 

predicted power-law dependence of the variance on the transmitter size as 

that for the incoherent source discussed in Section II.D. The reason for 

this is that in both cases, interference does not play a role: in the 

incoherent case, there is no interference, and in the coherent case, with 

W there is insignificant wavefront distortion. Realm (c) represents 

significant wavefront interference, a proliferation of diffraction scale 
26 

patches, and beamspread determined by diffraction from p Q . Realms (d,e) 
represent a saturated point source, and realm (f) represents a large 

source with multiple scattering. In the latter case, which has not been 
experimentally verified, the covariance scale is deduced from adjacent 
realms (c,e) by the requirement of consistency at their respective transi¬ 
tions. In all cases involving beam breakup or saturation, the normalized 
variance is unity in agreement with a model of random-phase, identical- 

3 

frequency summed oscillators. 


24. D.L.Fried and J.B.Seidman,"Laser Beam Scintillation in the Atmosphere", 
JOSA, 57_, 181, February 1967. 

25. A.Ishimaru, "Fluctuations of a Beam Wave Propagating Through p Locally 
Homogeneous Medium", Radio Science, 4_, 295, April 1969. 

26. J.R.Kerr and R.Eiss, "Transmitter-Size and Focus Effects on Scintilla¬ 
tions", JOSA, 62, 682, May 1972. 

27. J.R.Kerr and J.R.Dunphy, "Experimental Effects of Finite Transmitter 
Apertures on Scintillations", JOSA, 63^, 1, January 1973. 
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TABLE I. Focused TEM Q0 Beam (O 


/L/k< p 

r 




Table II: 

Collimated TEM^q beam. Realms (a,d,e) are identical to those in Table I. 
Realms (b,c) are essentially plane-wave illumination within the first- 
order theory, while (f) represents plane-wave illumination with saturation. 

Table III: 

This is the case analyzed in Sections II.A,B, in which a focused TEM Q0 

beam illuminates a diffuse target. The target is the "source" and the 

receiver is in the transmitter plane (Figure 1). "Near and far field" refer 

to the receiver-plane relationship to the illuminated spot on the target, 

which in some realms is in turn determined by the p diffraction scale. 

o 

Realms (a,b) are the first-order point source case with log normal scintil¬ 
lation statistics. In realms (d,e) the "atmospheric speckle" predominates, 
with joint gaussian field statistics, while in (c,f) the target speckle 
predominates, also with joint gaussian statistics. In realms (c-f) the 
model of independent additive oscillators applies, resulting in a normalized 

g 

variance of unity. In the transition region between (e,f), the field 
statistics are gaussian but not joint gaussian. In the transition region 
between (a,d), the normalized irradiance variance may exceed unity. The 
amplitude perturbation term (C^) is important (Section II.B) in realms 
(a,b). In fact, Figures 3a,c correspond to realms a (or b) and d, respec¬ 
tively. 

Table IV: 

Collimated TEM Q0 beam illuminating a diffuse target. This case is 
only of interest for >L/k, such that "collimation" has meaning (near 
field conditions). The basic considerations are the same as for Table III. 

(We also note here a further possibility, that of a direct monochromatic 
diffuse source (such as an extended laser transmitter with a diffusing 
screen in front of it.) In this case. Table I applies with the only 
exception that w = « and, in realm (b), o 2 = 1.) 

t n 

Table V: 

Incoherent (Nonmonochroinatic) Source. In this case we assume that 
the spectral width of the source is not sufficient for "spectral diversity 11 














III. Diffuse Target/Focused Beam (a ) 

o 

(Note: w = °°) 


vVk <p /L/k>p 

o o 



(a <p </L/k - Near Field) 
oo 


(a </L/k<p -Near Field) 
o o 




TABLE IV. Diffuse Target/Collimated Beam (a Q ) 
(Note; w - “) 


/L/k<p o /L/k>P Q 
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TABLE V. Incoherent Source (b) 
(Note: w = °°) 


b</L/k,p 

o 

(Far Field) 


/L/k<b<p 

o 

(Near Field) 

p <b</L/k 
(Far Field) 


b>/L/k,p 

o 

(Near Field) 




scintillations, i.e. smoothing owing to broadband effects. Unfortunately, 

the incoherent case is not yet perfectly understood, Realm (a) is a point 

source in the nonsaturated case. Realms (b,c) correspond to an extended 

source also in the single-scattering case, for which the analysis of Ref. 

20 predicts a -7/3 power-law behavior of the variance, and with a covariance 

scale that is related to the source size through geometrical effects (Ref. 

19 and Section II.D); since the source is incoherent, the condition p <b 

o 

in going from realm (b) to (c) has no effect. Realms (d,e) are for the 
saturated point-source, and it is expected that as the spectral width 
increases over that for the monochromatic case, the "atmospheric speckles" 
(P q ) are washed out, although the degree of spectral spread required to 
cause this is not understood. The variance will decrease from unity when 
the detector has insufficient bandwidth to detect the "beats" between the 
individual oscillators comprising the nonmonochromatic source. Further 
collapse is expected for the multiple-scattering, near-field case of realm 
(f). As mentioned previously, it would be desirable to explore (d-f) fur¬ 
ther using the "saturation covariance" work of Refs. 14rl7. 

II.F.Further Work and Applications 

The application of the efforts described in this report is to the 
performance of coherent optical adaptive (COAT) systems. The desirable 
analytical and experimentally verified results include the strength, spatial 
scale, and spatial and temporal spectrum of scintillations and their sta¬ 
tistics, for a variety of sources including single- and multi^-mode illumi¬ 
nated diffuse and structured targets. The coherent or complex amplitude 
correlation (mutual coherence function) of the receiver field is equally 
pertinent. 

Finally, the turbulence effects should be applied to control system 
models of actual adaptive systems which are generally phase correcting sys¬ 
tems. We have given this preliminary attention with the assistance of sys- 
28 

terns engineers, including the dynamics of target lock-on and the resultant 
evolution in spot and speckle conditions, 

Significant progress has been made on each of these problems, as 
described above. 

28. We acknowledge with appreciation the collaboration of James E, Pearson 
and others at the Hughes Research Laboratories, Malibu f California, 
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We detail below specific topics requiring further analytical work: 
Completion of Fundamental Understandings 

A. Complete the coherent-source, diffuse target covariance analysis 
of Section II.B., utilizing the full point-source mutual coherence function 
including C^ terms. This primarily reduces to thorough numerical analysis 
of the integral results; also, the "saturation" form of C (Refs. 14-17) 
should be incorporated. The results should be extended to the time-lapsed 
covariance, temporal and spatial spectra, and probability distribution of 
irradiance. 

B. Complete a similar analysis with glints and other target structure, 
i.e. through second-irradiance moments and spectra (Section II.C.), and 
including the complex amplitude (mutual coherence function). 

C. Incorporate incoherent receiver-aperture smoothing. 

D. Incorporate multimode sources, and the extreme of incoherent sources 
(Section II.D.) including the saturation form of C . I n particular, the 
incoherence or spectral width required to wash out^'atmospheric speckles" 

(Pq scintillation scale) should be understood. 

Applications 

A. Characterize target models at wavelengths of interest and incorporate 
them into the analysis. 

B. Analyze the effects of moving targets. 

C. Analyze the influence of all turbulence effects on the operation 
of real adaptive systems, including lock-on dynamics. 

Secondary Aspects 

A. Influence of < x <j>> or the phase-amplitude correlation in the point- 
source mutual coherence function. 

B. Effects of inner scale of turbulence, 

III. Experimental Effort 

The experimental effort is proceeding in three stages. The first 
involves preliminary experiments over a short (100m) path with a Helium- 
Neon laser operating at 6328 X, with a quasi-diffuse (Scotchlight) 
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This project was designed to permit shakedown of the experimental tech¬ 
niques and system electronics and verification of basic theoretical pre¬ 
dictions at an early date, and has been completed. The second ^cage 
involves operation over a 1 km uniform path with an Argon laser operating 
at 488C)R using a similar target. The required installations are nearly 
complete at our instrumented field site. The third stage involves the use 
of a TEA laser (3.5p and 10.6y) which is currently being fabricated and 
which will also be utilized at the field installation. Characterization 
of target models will also be conducted in the laboratory at wavelengths 
of interest, and target structure will be added at the field si*.e 

A basic problem in these experiments is the energy available off of a 
diffuse or quasi-diffuse target at a distance. At visible wavelengths, 
the directional properties of Scctchlight permit an adequate signal-to- 
noise ratio while preserving the pertinent properties of a "diffuse" 
target; that is, the speckles in the absence of turbulence have identical 
properties to those observed for a truly diffuse source. At the infrared 
wavelengths, the high instantaneous peak power of the TEA laser is required. 

The project goal is the observation and measurement of a number of 
statistical properties as described in the analytical discussion of Section 
II. The basic system is illustrated in Figure 5. The coherent laser 
illumination is modulated, focused, spatially filtered, and expanded to 
impinge on the remote target ~urface. The output beam may be focused, 
collimated, or operated in an intermediate condition. The return energy 
is detected by two separable receivers located adjacent to the transmitter, 
demodulated, filtered, and recorded on analog tape for digital processing. 
Simultaneously, the correlation of the two irradiance signals is computed 
electronically with variable time averaging and displayed. The measure¬ 
ments are repeated at various receiver separations to yield a covariance 
curve, and the turbulence strength is measured using a microtherma.1 prooe 
system. The output from final processing includes 

Turbulence level (C 2 and p ) 

n o 

Spatial covariance function 

Probability distribution of irradiance 


Modulator 
















Normalized variance of irradiance 
Temporal power spectrum 
Meteorological parameters 

As discussed in Section II, there are in general six possible parameter 
realms resulting from variations in transmitter configuration, wavelength, 
range, and turbulence strength. It is our intent to cover all of these 

realms; as many as possible will be covered at each of the operating wave¬ 
lengths . 

Initial Results 

The purpose for the short-path experiments was to aid in the design 
and shakedown of the optical, electronic, and data processing systems to 
be used in the further stages. The low path-integrated turbulence levels 
over such a path in the winter months result in large values of p compared 
to transmitter size (a q ) and Fresnel scale (L/k)* 5 . The relevant parameters 
are listed in Table VI, and it is noted that the basic parameter realm was 
limited to (L/k) ^ a Q < P Q . However, defocusing of the beam permitted the 
measurement of smaller speckles. 

Visual observations indicate the expected speckle pattern, and show 
that the pattern is nearly sationary. This indicates that the "target 
speckles" (a q ) are predominating. In order to achieve reasonable time 
averaging, it therefore proved necessary to rotate the target so that the 

speckle field is swept across the detectors. The rate of motion is limited 
by the receiver bandwidth. 

Covariance functions of irradiance have been obtained under a number 
of conditions. Figure 6a shows the measured covariance function for a 
focused beam over the 85m path length, with a = 0.013m, C 2 = 5.0xl0" 15 m 2/3 
and p Q = 0.15m. The theoretical curve is also shown, based on the reasonings 
of Section II.A. Since the point-source log amplitude variance (a 2 ) was 
quite low, with a comparable size for the target speckles and Fresnel scale, 
the amplitude perturbation term of Section II.B. was not pertinent (see 
Figure 3). Figure 6b shows a similar curve with a defocused beam such that 
the spot radius on the target is 5mm compared to 1.3mm in the preceding 
case. The expected decrease in speckle size is observed. 
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TABLE VI. Experimental Parameters 



Stage 1 

Stage 2 

Wavelength 

.6328p 

0.488p 

Range 

85m 

100-1000m 

Transmitter Beam Radius a 

o 

.013m 

.005-.02m 

/lTr 

.0029m 

.0028- 

•0088m 

/XL 

.0073m 

.007-.022m 

Range of values 

0.1*1.0 

meter 

1mm-10cm 
(including 
summer 
conditions) 

Laser Power 

^50 milliwatts 

^1.0 watt 

Modulation Frequency 

9 kHz 

100 kHz 

Signal Bandwidth 

1 kHz 

_> 1 kHz 

Receiver Size 

2 mm 

2 mm 

SNR (Scotchlight) 

o 

H 

A 

V 

H- ■ 

o 

-p- 

SNR (True Diffuse) 


10 3 


(Night con¬ 
ditions, 100m 
path) 
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Receiver Plane 


Target Plane 


"P i 
"?] 


Figure 7. Variables for mutual coherence function. 
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Long-time-averaged ('v 6 minutes) measurements of the normalized variable 
of irradiance are summarized below, for the focused case* 


C 2 (m _2/3 ) 
n 

a 2 /1 : 

- 1 ^ 

8.2x10 

1.12 

-1 ft 

8.9x10 

0,85 

4.5xl0 -15 

1.16 


The theoretical value is of course unity. 
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Appendix A. Generalized Mutual Coherence Functions 

The point-source mutual coherence functions are crucial to the utiliza¬ 
tion of the extended Huygens-Fresnel approach. In this Appendix we review 
the generalized, high-order mutual coherence function. The derivation will 
not be given but is based upon considerations suggested by Yura. 

The general n th order mutual coherence function is 

H(pi,P2 •••P2 n 5Pl »P2 ••• P2 n ) 

= <exp jtHPx,Pl)+ ^*(P2,P2)+---+ ’H^2 n _ r P'2 n _ 1 ) + ^*<P2 n .P2 n )( > (Al) 


where xp = X + 14* • 

We now make the following assumptions: 

1. x has a jo int normal distribution, with <x> = -® x ■ 

2. 4> has a joint normal distribution, with <4>> = 0. 

3. x £nd 4> are independent. 

There are limitations on the validity of these assumptions. It is known 
that x is not gaussian in the case of saturated scintillations, and as 
pointed out in the preceding report on this program, it is in general not 
jointly normal. 3 Also, < X 4» is in general nonzero and can be calculated 

from the first-order theory. 

Assumption #1 can be replaced with the assumption of a gaussian refrac¬ 
tive index field, bu<. this is a very poor representation of the actual 
nature of that field. 22 Alternatively, 2 we may simply say that the results 
to be given here are correct to "second order in the refractive index varia¬ 
tions", which is at best a heuristic approximation in conditions of saturated 
scintillations. However, we will not dwell further here on the degree of 
approximation involved in the assumptions, which comprises a fruitful area 
for further work. It is sufficient to point out that very useful results 

are obtained. 

The primary result is as follows: 
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J.R. Kerr, et al, " 
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Propagation of Multiwavelength Laser Radiation 
Turbulence", RADC-TR-74-320, November 1974 , 


(A003340). 
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where M is the number of j satisfying 2j+i=2n (for i even) and 2j+i=2n-l 

(for i odd). The wave structure function and log amplitude covariance 

function C comprise a generalization of those normally seen in the 
X 23 

literature. Referring to Figure 7, they refer to complex phase or log 
amplitude quantities as a function of (p^-p^) f- n the observation plane, 
for fields emanating from and p^. Specifically, 
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It may be noted that the functions satisfy (p,p) reciprocity. 

Special Cases 

(a) ._ n _ = l _ - Ij D 

H(Pi,P2;P;.P^) = e ^12 (A5) 


(b) n = 2 

H(Pi,P 2 ,P3,P 4 ;P1,P2,P a ,P4) 


e _l 5 (Dl2-D 13 +Di4+ D 2 3- D 24+ D 34)+2C^^+2 c ^^^ 
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This is a generalization (double-arguments) of the result of Ref. 23 


(c) 

(d) n 

H = e 


= 2: Pj = P 2 »P 3 = Pi»(calculation of <IjI 2 >) 

See Eqs. (18,20) in Section II. 

= 2; Pi=P2,Pa = P4;Pi=P2 = P3' = P4(calculation of <I^l 2 > 
source) 

’x (P3“Pl) 
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